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ABSTRACT 



A simplified form of the three-degrees-of -freedom nonlinear 
equations of the 3000 ton Surface Effect Ship (SES) was lin- 
earized around an equilibrium point when a thruster failure 
or thruster reversal occured. Optimal control feedback gains 
were calculated with the linearized system and applied to the 
nonlinear system. Evaluation of this design was based upon 
both comparison of the real-time human operator response and 
additionally a state feedback design based upon a trial and 
error procedure. 
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I . INTRODUCTION 



Over the past years a lot of attention has been given to 
the Surface Effect Ship (SES) due to its advantages over con- 
ventional ship designs [Refs, 1, 2, and 3] 

Recent investigations regarding water jet propulsion SES 
have been conducted on the computer assisted mode of operation 
either in a normal mode of operation or during a variety of 
failures [Ref. 3], Failures considered include, among others, 
complete loss of thrust and reversal of thrust of a nozzle. 

The purpose of this thesis was to investigate the possi- 
bility of optimal maneuvering of a 3000 ton Surface Effect 
Ship during the above mentioned emergencies . 

As discussed in [Ref. 3], turn maneuvers for the 3KSES 
are the result of turning moments generated by modulating the 
waterjet thrust vectors. This modulation is achieved by a 
combination of thrust amplitude modulation (controlled by 
varying the flow rates of the four water jets to create a 
net differential thrust, i.e., by varying the power turbine 
speeds of the individual propulsion engines) , and thrust 
vectoring (controlled by varying the directions of the four 
waterjet deflection sleeves). 

In the computer assist mode, the relative position of the 
helm determines the angular rate at which the turn is nego- 
tiated. Full helm position corresponds to the maximum safe 
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turning rate for a given ship speed. Calibration of the helm 
input is done within the control computer, which then gen- 
erates commands to the individual control effectors to produce 
the commanded turn rate . 

In the manual mode, the relative helm position generates 
commands for the direction nozzles directly. A common position 
command drives all four deflection sleeves. Because of that, 
the rate of turn which corresponds to a particular helm posi- 
tion is dependent on the hydrodynamic yawing moment which is 
a function of ship speed. The operator can also generate a 
thrust differential turning moment in the manual mode by 
proper selection of the propulsion engine control levers. 

To relieve the operator of the work load related to moni- 
toring and maintaining ship heading over extended periods of 
time, a heading control autopilot is included in the design 
of the heading control system. From an operational point of 
view, the heading autopilot is limited to a "heading hold" 
function, and is used only in conjunction with the computer 
assist mode of operation. The net result of this operation 
is that the heading autopilot replaces the operator's hand on 
the helm, while the helm is still the primary control input 
device for the control of ship heading. 

This is accomplished by direct feedback of the yaw rate 
to a processor which also gets an error channel from the auto- 
pilot through a control filter using proportional plus inte- 
gral control of the heading angle. 
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The fundamental optimality requirement used in [Ref. 2] 
is concerned with minimizing fuel and effector duty cycle 
while executing prescribed maneuvers. This was done by 
minimization of a function of the state variables. 

The approach used in this work was to achieve results by 

ir f T T 

minimizing a functional of the form / (x Qx + U Rp)dt 
to compute the gains for the state feedback control (see 
Appendix A) . 

The resulting linear regulator design was tested on the 
nonlinear system to compare the performance of the craft 
during a failure mode with the performance of a human operator 
using a real time simulation [Ref. 4]. The performance of the 
design during a failure was also compared with a corrective 
algorithm for control found by simple trial and error pro- 
cedures . 
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II. EQUATIONS OF MOTION OF THE THREE -DE'GREE- 
OF FREEDOM (3DOF) MODEL 



The three degrees of freedom to be used in this analysis 
are: surge, sway and yaw [Ref. 5]. Therefore, the following 

assumptions and approximations were made: 

1. Roll is zero. 

2. Pitch and draft remain constant. 

3. The origin of the ship's coordinate system is located 
at the center of gravity of the ship. 

4. At time = 0, the ship’s coordinate system and the 
fixed (geographical) coordinate system have the same origin 
and the same coordinate axes, i.e., x„.,, = x =0 and y. T/w , = 
y Q = 0 ( see Fig . 1) . 

5. At time = 0 the ship is heading in the positive x^y 
direction . 

6. Drift angle subtraction effect on thrust angle is zero. 

Then, as a good approximation, the equations of motion 

are : 



Surge : 


m(u-vr) = X 


(1) 


Sway : 


m(v+ur) = Y 


(2) 


Yaw: 


I r=N 


(3) 




z 


Navigation: 






x q = ucos^-vsin^ 


(4) 




y Q = usin^+vcosi^ 


(5) 
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X NAV 




Fig. 1 - Definition of Coordinate System 



NOTE: Rudder condition for right turn is shown. All runs 

start at geographic origin. 
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( 6 ) 



r 7 / 2 , 2 

V = / u +v 
s 

3 = tan“ 1 (-^-) (7) 

In the above equations: 

m = mass of the rigid ship 

I = moment of inertia about the z axis 
z 

u = V cos3 = velocity in the x direction (Surge) 
s 

v = V s sing = velocity in the y direction (Sway) 
r = ty = angular velocity about the z axis 
ip = heading (yaw) angle, Fig. 1 
3 = drift angle, Fig. 1 

5 = rudder angle, +6 produces a left turn 



Also : 




x = /x dt 
o o 


y Q = /yt 


u = /udt 


v = fv dt 


r = / rdt 


ip = /i|>dt 



The right-hand side of the equations of motion have been 
represented by X, Y, N. These symbols represent: summation 

of forces (X, Y) and moments (N) . Each summation consists of 
many terms. For the initial simplified model, only the few 
terms that are major contributions have been retained. Those 
terms are considered sufficient to produce a simulation which 
would be in close agreement with obtained data. At a later 
time, the model may be improved by adding more terms and ad- 
justing coefficients. 
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After reduction of the right-hand-sides, these force and 
moment summations become: 

X = GTcos6-C dx u 2 (8) 

Y = GTsin$-C DY v.|v.| (9) 

N = -(GTsin<$)l +C vlvll (10) 

o DY 1 1 w 



where (see Fig. 2): 

GT = gross thrust in lbs. 

(5=6 s = 6p = angle of all thrust vectors (starboard and port) 
in radians 

Cdx = coefficient of drag in the direction of surge 

C^y = coefficient of drag in the direction of sway 

1 = thrust level arm in feet 

o 

1 = sway-drag level arm in feet 



Also, note that the drag forces have been chosen as lumped 
parameters acting at the two points P S P', i.e. PSP* are 
equivalent force centroids. 

The length 1 is the distance from the center of gravity to 

the stern, and 1 is the distance (in the x-direction) from the 

center of gravity to the points PSP'. 1 is taken as a fixed 

parameter (it could be charged by various conditions of ship 

loading) but 1 is adjustable in order to be used as a para- 
w 

meter to make the model fit test data. 

The coefficient C^ x in eq. 8 is a lumped parameter equiva- 
lent for the hull friction effect (slender body theory). The 
2 

term C^u models the hull drag, aerodynamic and hydrodynamic 
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PORT 




Fig. 2 - Forces Acting For A Right Turn 



effects . The bubble drag is not included in the initial 
model . 

In like manner, the coefficient in eqs . (9) and (10) 

models the cross flow drag effects. Note also that the added 
mass term has been temporarily omitted and will be added 
later, especially in the 3KSES study since its contribution to 
the XR-3 model is almost nonexistent [Ref. 6]. 

Substitution of equations (8), (9) and (10) into the 
equations (1), (2) and (3) provides the basic equations of 
motion for the three-degrees-of -freedom model: 



1 

u = — 
m 


(GTcos<5 


• c dx u2) + vr 


(11) 


1 

v = — 
m 


(GTsin6 


- C Dy v j'v j ) - ur 


(12) 


n 

it 


■ (-(GTsin6)l + C nv v|v|l ) 
O JJI ’ 1 w 


(13) 



z 



Equations (11), (12) and (13) are applicable to any SES 
craft. Only the particular parameter values to be used in 
the equations would depend on the particular craft (XR-3 or 
3KSES ) . 

In the case of the 3KSES , the equations take the form 
( see Fig . 3 ) : 

u = - ( -C_. v u 2 + (T7+T8+T9+T10 )cos<5 ) + vr 
m DX 

v = - ( (T7+T8+T9+T10 ) sinS -C nu v ! v [ ) -ur 
m DY 1 1 

r = ~ (Cpy(ww)v | v | + [T7cos6 S^+T 8 cos 6 S^ 
z 

- T9cosSS^ - TIOCOS 6 S 2 - (T7sin6 +T 8 sin 6 + T9sin<$ 

+ T10sin6 ) (00 ) ] + A 90 uv(ww )) 
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= V 

x q = ucos^-vsinijj 
y Q = usin^ + vcosij; 

For the purpose of this study the thruster's angle 6 was 
the same for all thrusters except the starboard one, which 
in the case of thrust failure was considered 0° (thruster 
oriented straight ahead) and in the case of thrust reversal 
it was considered as being 180° (thruster oriented straight 
astern) . 

In the first case, the thrust of the failed thruster was 
set to zero pounds and in the second case the thrust was set 
to full value (98000 pounds for 56 Kts operation speed). That 
means that either in the case of thrust failure or thrust 
reversal of the starboard thruster, the rest of the thrusters 
are ganged together and act in the same direction. 

Therefore, the angle 6 can be redefined as Z7, Z8, Z9 and 
Z10 as shown in Fig. 3 to correspond to the appropriate thrus- 
ters and in both cases Z7-Z8-=Z9 with Z10 being either 0 or i 
radians . 

In order to linearize the equations, the following assump- 
tions were made : 

1. The maximum rudder angle is less than or equal to 30 
degrees and for most high speed manueuvers the rudder angle 
stays below 15 degrees. Thus, by letting cos6~l and sin6~6 
(where 6 is in radians), the difference from the actual values 
for 6= .2618 rads = 15° is 
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a . 


cosl5° = 


.966, difference = 


1 .966 6 x 100 


b. 


sinl5° = 


.2588, difference 


.2618-. 2588 
. 2588 



x 100% = 1.15% 

Differences with these magnitudes should not introduce 
large deviations from the linear design when applied to the 
nonlinear equations. Thus in the analysis <$ max was restricted 
to 15 degrees or if 6 x exceeds 15 degrees it probably will 
be of short duration and will then not introduce large errors 
in the design values for control. 

2. The yaw angle ip was considered to be less than 10 
degrees which agrees very well with the results obtained from 
the simulation of the craft. 

3. Even though the sign of the sway velocity changes de- 
pending on the particular situation occuring, the term v[v| 
was considered as v 2 . 

With the above assumptions, the equations for the 3KSES 
become : 

f, = u = - ( -C tw u 2 +T7+T8+T9+T10 ) +vr 

1 in DX 

f 0 = V t - ( (T7+T8+T9+T10)6-C nv v 2 ) - ur 

2 m DY 

f 3 = r t (C dy v 2 (ww)+T7S 4 +T8S 3 -T9S 1 
z 

- T10S 2 -(T7+T8+T9+T10)6 (00 )+A 22 uv(ww) ) 

= 'P = r 

f 5 = *o ^ U " V ^ 
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f 6 = y 0 : uy + v 



Defining the state variables 



x, =u-u 
1 o 



x 2 = v - v o 



x 0 -r-r 
3 o 






X =X -X 

5 o oo 



x =y -y 
6 J o oo 



y, =T7 -T7 
1 o 



y 0 =Z7-Z7 
2 ( 



Then, x^ = u 



x 2 = v 



x 3 = r 



x 4 = * 



Xr = X 

5 o 



x = y 
6 J o 



The linear state equation coefficients are 



3f 

lu" 



1 _ 



■ 2C dx u 



m 



u=u 



- 2C PX u o 

m 
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3f 



1 



3v 



r 



r=r 

o 



v 

o 



3f ] 

3r 



V 



v 



v=v. 



3f 1 af x afj 
a "ax ay 



= o 



a f i 

3T7 



1 

m 



3f i 

az7 



= o 



af, 

lu" 



= -r 



= -r 



r=r 



a f ; 

"av" 



- 2C dy v 

m 



v=v 



~ 2 C DY v o 

m 



af. 

W 



= -u 



= -u 



u=u 



af„ af„ af. 



aip 



3x ay 

o J o 



= 0 



3f 2 

"3T7 



Z_7 

m 



Z7 = Z7 



Z 7 



m 



23 



3f 



2 T7+T8+T9 



3Z7 



m 



T7 = T7 C 

T8=T8 

T9=T9 



T7 +T8+T9 
o 

m 



9 f 3 A^vCww) 



3u 



A 22 V o ( ww) 



v=v 



3f , 
3v 



2 C D y v ( ww ) + A 2 „ u ( ww ) 



<2C DyV A 22 U o )( ww) 



u = u 
v=v c 



3f 



3 _ 



9f 3 . 9f 3 . 9f 3 



3r 



3^ 3x 3y 
o o 



9f 3 

3T7 



S 4 -Z7(00) 



S.-Z7 (00) 
4 o 



Z7 = Z7 



3f 



3 -(00) (Z7+Z8+Z9) 



3Z7 



- ( 00 ) ( Z7 +Z8+Z9) 
o 



Z7 = Z7 



3 f 4 3 f 4 



af 4 9f 4 3f 4 3f 4 



3u v 3 <Jj 



3x ~ 3y 3T7 3Z7 

o J o 



= 0 



3f L 

3r 



= 1 



3f r 



3u 



24 



3f 



5 



3v 







-Ip 



o 



3f 



5 



dip 



-v 



v=v 

o 



-V 

o 



3f r 3fr 3f, 3f c 3f c 

5 _ 5 _ 5 _ 5 _ 5 

3x 3r 3y 3T7 3Z7 

o J o 



3f 6 

3u 




* 



o 




af„ 

6 



3f _ 3f _ 3f„ 3f _ 3f . 

6 6 _ 6 _ 6 _ 6 _ n 

3r ~ 3x (V ~ 3T7 3Z7" U 

o J o 

In standard form the system becomes: 
x = Ax + B y where 

r* X ~ 



u 



u = u 
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* to 



X 



1 



« 

X 



2 



% 

X 




X 



5 



* 



6 



u 



AT 7 
AZ 7_ 



a ll 


a 12 


a 13 


0 


0 


0 


1 — 1 
CM 

rtf 


a 2 2 


a 23 


0 


0 


0 


a 31 


a 3 2 


0 


0 


0 


0 


0 


0 


a 43 


0 


0 


0 


a 51 


a 5 2 


0 


a 54 


0 


0 


l — 1 

CD 

rtf 


a 6 2 


0 


a 64 


0 


0 



'll 



'21 



'31 



0 

0 

0 



22 



32 



0 

0 

0 
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and 



l ll 



2C PX U o 

m 



a 12 = r c 



13 



v 



a 21 = “ r c 



22 



2C DY V o 

m 



a 23 = " u c 



31 



3 2 



A 00 v (ww) 
1 1 o 



(2C DY V o +A 22 U 



a 43 ” 1 



51 



a 5 2 = “*< 



a 54 = “ V c 



a 61 = K 



a 6 2 = 1 



a c „-= u 
64 c 



) ( ww) 
o 



27 



'll 



1 

m 



Z7 

c 

*21 " m 



22 



T7 +T8+T9 

o 

m 



31 



32 



S.-Z7 (00) 
4 O 



- ( 00) (T7 +T8+T9 ) 
o 



So, by linearizing around the equilibrium point (u q , v q , 

r , ip , x , y ) the equations become: 
o r o oo J oo 

-ZC-.yU , 

x, = — x, + r x 0 + v x_ + — C.. 

1 m 1 o2 o3ml 

2C nv v 27 

& 9 = -r r, X 1 “ ~ ~ X 9 " u ~ x o + ~~Zr~ 

2 ol m 2 • o 3 ml 



T7 +T8+T9 



+ (- 



m 



-) C, 



x. 



A 22 v o ( ww) ^ <2 C dyV A 22 u o H wv,) 

1 X 1 + 1 x 2 



S, -Z7 (00) ( 00) (T7 +T8+T9 ) 

+ ( 4 ° ) c 5 c 

k I J C 1 I 2 



X 4 = X 3 



x c = X.. - ip x 0 - v x,, 
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The above equations were used to determine by computer 
simulation the required gains for the state feedbacks using 
the linear regulator solution [Refs. 7, 8 and 9]. 

In order to obtain the equilibrium point, the nonlinear 
equations were used by setting u=0, v=0, r=0 and determining, 
after a failure had occured and the system being corrected as 
far as was possible by the trial and error algorithm (see 
Chapter III), the ''steady state" values for the above men- 
tioned variables. 

Two "steady state" runs were made, one with the starboard 
thruster failed and the other with the starboard thruster 
reversed . 

The "steady state" values obtained are listed in Tables 
I and II. 



TABLE I 

Steady State Values (Thruster Failure) 

Surge Velocity Sway Velocity Angular Velocity Yaw Angle Track 

Deviation 



( ft/s ) 


(ft/s) 


(rad/s) 


(rad) 


(ft) 


74.422 
Port Thrust 
(lbf ) 

87600 


1.86 


-. 7415xl0~ 7 


-.0254 


10.88 
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TABLE II 



Steady State Values (Thruster Reversal) 

Surge Velocity Sway Velocity Angular Velocity Yaw Angle Track 

Deviation 

(ft/s) (ft/s) (rad/s) (rad) (ft) 

55.66 2.62 -.3417xl0~ 5 -.04709 21.41 



For this study, the steady state value achieved by the 
simulation with the thrust failure were used for the develop- 
ment of the linear system and the subsequent use of the opti- 
mal control scheme. The choice was quite arbitrary but, 
nevertheless , the results show that the system does behave 
well, no matter which one of the two emergencies (thrust 
failure or thrust reversal) occurs. 
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III. TRIAL AND ERROR DESIGN 



A. INTRODUCTION 

Taking into account the impressive speed of the surface 
effect ship, it is evident that there is a need for controlling 
it either during straight runs or during an emergency which 
would be manifested by the loss of thrust of one thruster or 
even by the thrust reversal . 

As was seen from the results of computer runs with the 
model, the behavior of the ship during a failure was pre- 
dictable but nevertheless disturbing because the ship veered 
quite rapidly from the track that it was following before 
the failure. 

Therefore, the goal of the control procedure was twofold. 
First, it was required to maintain control of the ship and 
second to maintain that control in some acceptable manner. 

B. ANALYSIS 

As was evidenced (see Fig. 4), when the outer starboard 
thruster experienced a complete thruster failure, the ship, 
without any corrective action with respect to nozzle movement 
or thrust variation of the remaining jet nozzles, veered to 
the right and actually entered a turn. 

By analyzing this behavior, it was seen that the para- 
meters changing the most were the yaw angle ip and the deflec- 
tion from the original track which was assumed straight ahead. 
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Intuitively then, it was evident that from the available 
variables, ip and y Q would be most eligible to control in order 
to bring the ship back or close to the original track. Be- 
cause it was quite clear that, if the heading (yaw angle) was 
minimized as well as the deviation from the track, the ship 
would follow again the original track or a track that was 
parallel to that one. 

Therefore, the corrective action chosen was nozzle de- 
flection as a function of the yaw angle and the deviation 
from the original track. 

Furthermore, since positive nozzle angle drives the ship 
to the left, by selecting deviation to the right as being 
positive and yaw angle to the right being positive to the 
right, it was quite clear that the nozzle action would only 
involve ip and y Q without any sign change. 

Thus , the problem reduced to the point of determining the 
appropriate gains for ip and y^ and the rudder action took the 
form z=K-^^+K 2 y Q . 

C. TRIAL AND ERROR PROCEDURE 

As a beginning, was arbitrarily selected to be 1000 
and ]<2 being 5 , It was seen that the action of the rudder 
was too heavy and the movement of the ship was quite erratic 
and at the end became unbounded . 

Therefore, both gains were decreased to K^=8 and j< 2 =. 035 , 
which gave a more acceptable performance (Fig, 14), 
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The thus established rudder algorithm proved to be adequate 
even when a reversal failure of the starboard thruster was 
inflicted (Fig. 21). 

As a further improvement, the idea of modulating the 
thrust of the left thruster during a failure was brought into 
effect. That is, the thrust was varied to smaller values by 
the amount of the deflection of the ship off the original 
track ( Fig . 20 ) . 

The above scheme was used for all the runs presented in 
Figs. 14 through 27. 

In order to evaluate all the above corrective actions , 
the performance indices described in Chapter IV were used. 

As can be seen, the ship actually moves down the track 
with a large drift angle. That is, although a corrective 
rudder action is maintained, the ship does not turn into a 
circle and also it does not veer to the left or to the right 
although there exists a sway velocity and asymmetrical thrust. 
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IV. LINEAR REGULATOR DESIGN 



A. INTRODUCTION 

The theory behind the linear regulator is well established 
and need not be repeated here. For a brief discussion, see 
Appendix A. 

After the linearization of the system's equations around 
an equilibrium point, the IBM/ 3 60 was used to obtain the 
solution of the matrix Riccati equation, using the method 
presented in [Refs. 8 and 9]. 

B. ANALYSIS 

It was evident that in order to achieve a parallel track 
for the craft after a failure had occured , the deflection 
from the original track should come into play. Also, in 
order for the craft to maintain this parallel track, the yaw 
angle should be minimized as well as the angular velocity 
about the z axis . 

As with the trial and error procedure the corrective 
action chosen was ganged deflection of the operating nozzles 
and thrust control of the port thruster to compensate for the 
failed starboard thruster. 

C. PERFORMANCE INDICES 

One of the basic problems of modern control theory is 
the translation of system specifications, often in such 
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subjective terms as "good rise time with reasonable overshoot", 
into the form of a performance index. Furthermore, after the 
selection of a performance index, the designer has to make 
an initial guess for the weighting of the state variables and 
adjust them after comparing time histories of the controlled 
system . 

For the purpose of this study, the following performance 
index was selected: 

J = / t f (x T Qx + y T Rp) dt 
o 



Concerning the weighting of the variables, considerable ex- 
perimentation was done. In the final form of the above 
performance index the weighting matrices Q and R were: 



Q = 



0 

0 



0 

0 



0 0 
0 0 



0 

0 



0 

0 

20000 

0 

0 

0 



0 

0 

0 

8 

0 



0 

0 

0 

0 

0 



0 0 



0 

0 

0 

0 

0 

.035 



R = 



1 

0 



0 



1 
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The final weights used reflect the importance given to the 
yaw angle, the angular velocity about the 2 axis and the 
craft deviation from the original track. 

The optimal control inputs obtained are: 



0 0 



y 






-.0316 



.684 



-.00046 -.0007 0 

185.7 110.7 0 



0 

.186 



x 



The response of the nonlinear system with the above 
feedback is shown in Figs. 28 through 41. 

In order to compare the performance of the craft under 
the trial and error algorithm, the optimal algorithm and the 
real time simulation, the following two additional perfor- 
mance indices were used: 



J 



1 



f^f y dt and J 0 
o 2 

o 



/ t f ( Z 7 ) 2 dt . 
o 
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V. COMPARISON OF RESULTS 



A. THRUST FAILURE MODE (OPTIMAL ALGORITHM VS. TRIAL 

AND ERROR ALGORITHM) 

Observation of the sea track (x vs . y ) (Figs. 14 and 28) 
indicates that the trial and error algorithm is much less of 
the desired control algorithm when compared to the optimal, 
because it appears that the trial and error algorithm does 
not react fast enough to stabilize the vehicle whereas the 
optimal does so immediately. The trial and error algorithm 
does seek out the steady state value of the track deflection 
but the comparison of the two tracks is dramatic in showing 
how much better off is the optimal algorithm in controlling 
the craft during a thrust failure. It must be kept in 
mind in making this observation that the trial and error 
design used in this comparison was not developed to its full 
capabilities. As stated earlier, the trial and error design 
was intended only as a means of establishing the linear 
system design. 

From the time history shown in Figs. 14 and 28, the vehicle 
for the same failure scenario reaches a steady state point with 
the optimal control in 22 seconds whereas with the trial and 
error control the steady state point would be reached some- 
where around 200 seconds or even greater, due to the lightly 
damped control system. 
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As can be seen from Figs. 20 and 35 the port thruster is 
much more active with the trial and error algorithm. The 
range of modulation covers about 14000 pounds whereas utiliza- 
tion of the optimal algorithm gives a total thrust modulation 
of 2 pounds. This does not reflect exactly what happens to 
the port thruster, because as stated earlier in Chapter II as 
soon as the failure occurs, the port thrust drops to 90% of 
its full value . 

It should also be noted that the trial and error solution 
seeks the steady state value of the port thruster obtained by 
the optimal solution. 

Comparing Figs. 18 and 32 it can be seen in Fig. 32 that 
the forward and side thrust as well as the yaw moment reach 
a steady state within 10 seconds after the initiation of the 
thruster failure when the optimal solution is used. Instead, 
when the trial and error solution is used, the craft as shown 
in Fig. 18 seeks the same steady state but even at 100 sec 
after the initiation of the failure it has not stabilized. On 
the average, the values obtained from both algorithms are 
approximately the same. 

However, it is readily apparent that the trial and error 
control system has poles close to the joj-axis and as a result 
exhibits a response of a lightly damped system whereas the 
optimal control system responds with much faster time constants, 
which would indicate that this system has poles at a greater 
distance from the jw-axis. 
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B. THRUST REVERSAL MODE (OPTIMAL ALGORITHM VS. TRIAL 

AND ERROR ALGORITHM) 

All the previous comments that were made for the thrust 
failure mode also apply to the case of a thrust reversal for 
the starboard thruster. However, due to the increased 
severity of the forces acting during the thrust reversal and the 
greater magnitude of the yaw moment (see Figs. 32 and 39), 
the optimal control exhibits a slightly longer settling time 
(Figs . 28 and 35) . 

During the thrust reversal scenario, the absolute magni- 
tudes of the state variables is larger than those experienced 
during the thrust failure case. 

Comparison between the optimal solution time histories 
and trial and error solution time histories again confirms 
the fact that optimal state variable feedback provides a more 
highly damped system, with steady state values being reached 
within 20 seconds from the initiation of the failure (see 
Figs. 36, 38, 39 and 40). The trial and error solution is 
shown in Figs. 22, 24, 25 and 26 which again exhibits a 
lightly damped time history which would not be acceptable 
either from a mechanical or ship's efficiency viewpoint. 

It is observed (Fig. 23) that the trial and error solution 
requires the operational thrusters to be extended to their 
mechanical limits of +.5 radians even at the time =100 seconds 
point, whereas the optimal solution in Fig. 37 has a peak 
value for the thrusters of approximately +.4 radians and 
never requires them to extend to their mechanical limit. 
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Fig, 4 - Sea Track, Thruster Failure, No Correction 



40 



xlO (ft) 




(oas/q.j)u 

OS *5 QC ‘C 09 ‘0- 02 ' I - 03 ‘ I- 0 t‘i- GO'S- 09 •£- 02 *fi- 

(09S/4.J)A 

Fig. 5 - Surge and Sway Velocities Vs. Time, Thruster Failure, 
No Correction 



41 



TIME ( sec ) 



ca U 




CS *' 1 Oft * I 02 “I 00 * i 0S*0 CS*G Oft ’0 02*0 00 *C 

(spea)(fi 



‘S ‘0 9S 'C Sn "C OfTC c£ *C fi2 * C £ 

(spea) oixg 



30-0 GC ' 3 



Fig. 6 - Yaw Rate, Heading and Slip Angle Vs. Time, Thruster 
Failure, No Correction 



42 



10. on 30.00 30.00 MO. 00 50.00 00.00 /(). 00 80.00 90.00 100.00 

TIME(sec) 



a) 
a 
U 
o 

a) 

TJ 

cn T 



0) 

a 

sp 

o 

Ph 

SP 

(U 

s 

IP 

o 

Ph 



0$ 



O S 3 



© S3 



+J 

C 

a) 

S 

o 

s 

c 

IP 

3 



o 

o 



_o 

<x 



o 

o 

_o 

ao 



o 

o 



a 

h Sw 



o 

o 



o 

o 

l-° 



o 

o 



c 

o 



ocoos oo*GSf' oo • ooti oo - ose oo’ooe oo - os2 oo - oo2 

(JQD^OIX 30 ^ 0 ! pa^waoj 



00 ‘0S i oo - ooP 



C2 - ! 05 " 0 0S‘0 0t’0 00‘0 0£*0- 09’0- 

(jqi) soaoj 9 pT S 



os’o- 



02 ‘ I - 



CC '093 00 ' 008 OO'Cci 00 * 0Ti9 00‘GSS 00 * 08Ti OC'OOfi 00'02£ 0C‘0ri2 

(qj‘ jqi) , OI^quouioH uani 

n 

Fig. 7 - Forward Force, Side Force and Turn Moment Vs. Time, 
Thruster Failure, No Correction 



43 




si *0 ocrc- sro- oe*c- sti *o- 09‘c- sz. - o- 

( OSS/^J)A 



OS'O- 



SO ‘ I - 



OS *0 



Fig . 8 



OS’O Ofi’O OS ‘ 0 02 'C 01 '0 OC'O- 

( ^sas/spea) -,_oixa 

- Surge, Sway and Yaw Acceleration Vs. 
Failure, No Correction 



ox * c - 02*0- 

Time, Thruster 



44 




■s s a-. 



UJ) T 0IX x 

Fig. 9 - Sea Track, Thruster Reversal, No Correction 



45 



.00 MO. 00 80.00 120.00 160.00 200.00 2M0.00 280.00 320.00 360.00 MOO. 00 

y xl0 1 (ft) 




08 ‘o' 



00*0 



08*0- 



09*1- 



0ft*S- 00*8- 0 0 * ti - 

(oas/q-j) a 



08 ‘"ft- 09*5- 



Fig. 10 



Surge and Sway Velocities Vs. Time, Thruster 
Reversal, No Correction 
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Forward Force, Side Force and Turn Moment Vs, Time, 
Thruster Reversal, No Correction 
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No Correction 
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Surge and Sway Velocities Vs, Time, Thruster 
Failure, Trial and Error Solution 
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Fig. 16 - Effector Action Vs. Time, Thruster Failure, Trial 
and Error Solution 
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Fig, 20 - Port Thruster Action vs. Time, Thruster Failure, 
Trial and Error Solution 
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Solution 



57 



IJ.00 8.00 12.00 16.00 20.00 2*1.00 28.00 32.00 36.00 *40.00 



00 'S6 



00 '06 




0S‘£ 00 ‘ 9 OS ‘ ii 00 ‘ £ OS ‘ I OO'O CS't- 

(D3S/aj)A 



00 ‘£- 



OS’fi- 



Fig. 22 



Surge and Sway Velocities Vs. Time, Thruster 
Reversal, Trial and Error Solution 
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Fig. 23 - Effector Action Vs. Time, Thruster Reversal, Trial 
and Error Solution 
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Yaw Rate, Heading and Slip Angle Vs. Time, Thruster 
Reversal, Trial and Error Solution 
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Fig. 26 - Surge, Sway and Yaw Accelerations Vs. Time, Thruster 
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Fig. 27 - Port Thruster Action Vs. Time, Thruster Reversal, 
Trial and Error Solution 
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Fig. 28 - Sea Track, Thruster Failure, Optimal Solution 
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Surge and Sway Velocities Vs. Time, Thruster 
Failure, Optimal Solution 
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Fig , 30 - Effector Action Vs, Time, Thruster Failure, Optimal 
Solution 
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Fig. 34 - Port Thruster Action Vs. Time, Thruster Failure, 
Optimal Solution 
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Surge and Sway Velocities Vs, Time, Thruster Reversal, 
Optimal Solution 
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Optimal Solution 



73 



0.00 10.00 20.00 30.00 ‘10.00 50.00 60.00 70.00 80.00 90.00 100.00 

TIME ( sec ) 




(oas/ppa) ^._0 IXlX 



91 ‘0 80 "6 00 ’C- 8C *0- 91 ‘0- fi2*0- 28 ‘ 0 - 0 ti * o' - 8 fi * o - 

(P pa ) I _OT x 'f l 

91 '0 80 "0 00 *c- S0 # 0- SI * 0 - ti2 *c- S£ "0- OVC- 8Ti *0- 

(P*0C) 0TX8 



Fig. 38 - Yaw Rate, Heading and Slip Angle Vs. Time, Thruster 
Reversal, Optimal Solution 
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The sea track (see Fig. 35) of the optimal solution 
becomes stabilized at approximately 18 feet to the right of 
the original track (due north) , with a slight tendency to 
drift to the right. This is not of significant consideration 
due to the relatively small magnitude of drift and the fact 
that during the conduct of an actual emergency situation 
proper corrective action would be also taken by the crew of 
the ship in the first 10 to 15 seconds of the emergency. 

It is important to keep in mind that the ship's response 
for both the trial and error and the optimal solutions was 
the result of totally automatic operation and any human 
operator inputs were not considered. 

The deviation to the right obtained by the use of the trial 
and error algorithm at 100 seconds as shown in Fig. 21 has an 
average value of about 18 feet and is fluctuating around this 
value by a magnitude of 6 feet which would be more than enough 
to discrupt normal ship activity. 

In addition, the situation is more aggravated when surge, 
sway and yaw accelerations (see Figs. 26 and 33) are taken 
into account because even though they are small , they are 
still significant after 100 seconds. 
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C. COMPARISON OF THE OPTIMAL SOLUTION AND THE 

REAL=TIME HUMAN OPERATOR RESPONSE 

As stated in the introduction, a real time simulation 
of f ive-degree-of-freedom Surface Effect Ship was available 
[Ref. 4], utilizing a human operator to input only ganged 
thruster action, the thrust being constant. 

Comparison of the optimal solution and the above men- 
tioned simulation gave the following results during a thrust 
reversal scenario: 

1. Initial recognition of the emergency by the human 
operator did not occur until approximately three seconds into 
the emergency, as .compared to the optimal solution with an 
instantaneous recognition of the failure. Consequently, 
as shown in Table IV, the absolute values of all the outputs 
reached by a human operator controlled SES are much greater 
than those observed in the optimal control, given the same 
initial conditions. 



TABLE III 

Comparison of Optimal Solution and 
Human Operator Response 

HUMAN OPERATOR OPTIMAL SOLUTION 



Time for first overshoot ( sec) 


39 


22 


Deflection from original 
track (ft) 


100 

(worst case) 


18 


Time to settle back(sec) 


78 


35 


Final deflectionCft ) 


0 


18 


Turn rate (rad/sec ) 


.007 


.0065 
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Specific examples are as follows: 

1. Maximum track deviation obtained by the ship with the 
optimal control in effect was approximately 18 feet, whereas 
with the human operator maximum track deviation varied from 
35 feet to 100 feet, depending upon operator response and 
familiarization with the failure as well as the ability of the 
operator to recognize the onset of an emergency situation 
through the visual cues provided by the real time simulation. 

2. Turn rate was approximately the same, i.e., .007 rads/ 
sec for both situations. However, as might be expected from 

a human, the turn rate exhibits a much more sinusoidal be- 
havior as compared to the optimal solution which immediately 
sought and settled down to a steady state value. 

3. The state variables under real time simulation showed 
a more pronounced fluctuation. 

4. An advantage of the real time simulation over the 
optimal solution was that the human operator in all cases, 
i.e. both thrust failure and thrust reversal, was able to 
maneuver the ship back to the original track, whereas the 
optimal solution maneuvered the ship to a parallel track and 
never returned to the original. This, of course, stems from 
the method used to find the "steady state" values around 
which the equations of the system were linearized. 

Nevertheless, the difference of 18 feet between the original 
track and the one that is reached at steady state by the opti- 
mal algorithm is not significant compared either to the ship's 
width (100 feet) or to the tactical considerations of the ship. 
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D. NONLINEAR VERSUS LINEAR SYSTEM RESULTS 



In order to determine how "optimal" the response of the 
nonlinear system was with the use of the state feedback 
values obtained from the linear regulator design, a com- 
parison of values of the states and the performance indices 
were made during a failure. 

The results are shown in Table III. As can be easily 
verified, the differences in steady state values obtained 
from the linearized system and the nonlinear system are 
quite small, ranging from .77% to 1.77%, with the exception 
of the steady state thruster angle (15% difference) . 

The above results point out that the linear assumptions 
and approximations used in this study hold quite well over 
the operating range. 



TABLE IV 



Comparison of Linear and Nonlinear Models 



LINEAR NONLINEAR 



Surge velocity (ft/s) 
Sway velocity (ft/s) 
Angular velocity (rad/s) 



75.98 74.63 



L. 8 8 1.87 

. 214xl0"" 5 -.211x10 



Yaw angle (rad) 



-.0248 -.025 



Displacement in x direction (ft) 
Displacement in y direction (ft) 
Steady state thruster angle (rad) 



7453 7342 



10.762 10.88 



,15256 .1804 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



A, CONCLUSIONS 

Comparison of the results of all three methods examined 
leads to the conclusion that the optimal control solution is 
by far the best one. 

Control effort, be it thruster angle or thruster magni- 
tude variation, is less with the optimal control solution. 

This is of great consideration not only from the aspect of 
wear and tear of the mechanical parts but also from the 
point of ship performance as a naval unit, especially since 
its speed is so far greater than the conventional ships . 

Augmentation of human operator control inputs with the 
optimal solution during failure scenarios presents the best 
solution in order to maintain positive control of the ship 
under any operating conditions. This approach would combine 
the best attributes of the human operator and the optimal 
control schemes. 

B. RECOMMENDATIONS 

Based on the previous analysis , the following could be 
investigated in the future. 

1. A sufficiently complex constrained parameter optimi- 
zation routine (COPES-CONMIN) can be used to arrive logically 
at a set of relative system response with respect to physical 
constraints. This optimization will still be dependent on the 
performance index selected by the design engineer and his 
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his ability to calculate explicitly its value for a given 
dynamic response [Ref. 11]. 

2. The right-hand side of the equations of motion could 
be expanded to include more terms to provide an increase in 
quantitative accuracy of the drag terms such as bubble drag 
and include also a constant power propulsion plant versus a 
constant thrust one. This could be done easily since the 
terms are additive. 

3. Development of a computer program to augment the 
control efforts commanded by a human operator should be 
investigated, tested and implemented on the RTS5D of Ref. 4. 

4 . It might prove worthwhile to investigate more closely 
the trial and error design for improved response character- 
istics of this method of control. 



APPENDIX A 



LINEAR REGULATOR DESIGN THEORY 

Optimal control theory provides a means for utilization 
of state variable feedback in order to control a system by 
minimizing a performance index or cost function [Refs. 7 and 
8 ], 

Having a linear system defined as x = Ax + By where y 

is the control input vector, the linear op imal control law 

is defined as y , (t) = -Gx(t). 

-opt 

The assumptions under which this solution is valid are: 

1. The control vector is unconstrained. 

2. The performance index is defined as 

J = /"*■ f (x^Qx + y^Ry) dt 
o 

where Q is a positive semidefinite or positive definite real 
symmetric matrix and R is a positive definite real symmetric 
matrix. The matrices Q and R determine the relative impor- 
tance of the square of the error and energy expenditure 
during the transient period. 

3. The system is controllable and state observable. 
With these assumptions, substitution for y leads to 

J = / ^ (x^Qx + G ( t ) T x T R G(t)x) dt 
o 

J = / t f{x T (Q + G ( t ) 1 R G(t) )x> dt 
o 
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The matrix G(t) is found to be 

-v 

G(t) = R _1 B T P(t) 

*+ rw 

Here P(t) satisfies the matrix Riccati equation 
P(t) + Q - P(t) BR -1 B T P(t) + P(t)A + A T P(t) = 0 

subject to the boundary condition P(tf) = 0. 

The elements of the matrix G(t) are referred to as the 
feedback coefficients, since the optimal control consists of 
a time-weighted linear combination of the state variables. 

This result is important since it indicates that in an 
optimal system all the state variables are to be fed back, 
not only the output, as it is customary in classical control 
theory. A graphical representation of this optimal system 
is presented in Fig. 42. 

The Riccati equation is integrated backwards in time from 
the known terminal condition P( °° ) = 0 until the desired 
degree of accuracy is reached. When the final time approaches 
infinity, the P matrix becomes constant, thereby simplifying 
its determination. 

As can be seen for the previous presentation, a different 
linear optimal control is generated for every relative set 
of weighting matrices, Q and R, provided. Additionally, the 
system response is optimal in a mathematical sense only, since 
no assurance on acceptable system overshoots and input peak 
values can be guaranteed. 

A minimum area could have an excessive peak which would 
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provide a physically unrealizable control scheme. The 
design synthesis thus requires an initial guess at the rela- 
tive importance multipliers, followed by subsequent adjust- 
ments based on time simulation of the controlled system 
dynamics . 




Fig. 42 
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